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METHODS AND APPARATUS FOR TRUNCATION 
COMPENSATION 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. provisional 
application number 60/416,072 filed October 4, 2002, which is hereby incorporated in 
its entirety. 

BACKGROUND OF THE INVENTION 

[0002] This invention relates generally to methods and apparatus for 
computed tomographic (CT) image reconstruction, and more particularly to methods 
and apparatus that provide for truncation compensation. 

[0003] Under some scanning conditions, portions of a patient may 
extend beyond a region measured by a detector, which may lead to image artifacts and 
an incomplete representation of the imaged object. Some known methods have been 
published that address artifact reductions but not the imaging of the portion of the 
patient that is outside the field of view (FOV). However, it is desirable to image the 
portion of the patient that extends beyond the FOV such as is described in copending 
application serial number ##/###,### titled Methods and Apparatus for Truncation 
Compensation, and filed on June 24, 2003 with attorney docket number 129993 
(which is also hereby incorporated in its entirety), and claiming the benefit of 
provisional application serial number 60/416,072 filed October 4, 2002. Additionally, 
it is desirable to increase image quality of portions of the patient that extend beyond 
the FOV. 

BRIEF DESCRIPTION OF THE INVENTION 

[0004] In one aspect, a method for determining whether a projection 
is truncated is provided. The method includes calculating a sum of all samples at each 
projection view of a scan of an object, determining a maximum value of the calculated 
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sums, averaging a plurality of samples m at a projection view index k when the sum of 
all samples at index k is less than a predetermined percentage of the maximum value, 
comparing the average to a threshold t, determining the projection truncated when the 
average is greater than t, and determining the projection not truncated when the 
average is not greater than t 



imaging apparatus includes a radiation source, a detector responsive to radiation 
positioned to receive radiation emitted from the source, and a computer operationally 
coupled to the radiation source and the detector. The computer is configured to 
calculate a sum of all samples at each projection view of a scan of an object, 
determine a maximum value of the calculated sums, average a plurality of samples m 
at a projection view index k when the sum of all samples at index k is less than a 
predetermined percentage of the maximum value, compare the average to a threshold 
t, determine the projection truncated when the average is greater than t, and determine 
the projection not truncated when the average is not greater than t. 



which is encoded with a program configured to instruct a computer to calculate a sum 
of all samples at each projection view of a scan of an object, determine a maximum 
value of the calculated sums, average a plurality of samples m at a projection view 
index k when the sum of all samples at index k is less than a predetermined 
percentage of the maximum value, compare the average to a threshold t, determine the 
projection truncated when the average is greater than t, determine the projection not 
truncated when the average is not greater than t, and estimate a total attenuation x(k) 

kn — k k — k\ 

in accordance with r{k) =— £(£i) + — £(^2) where ki and k 2 are averages 

of a plurality of view locations of un-truncated views adjacent to a truncation region 
comprising at least one projection determined truncated, and £(k) is calculated as 



7 P(Uk) • The program is further configured to instruct the computer to 



[0005] In another aspect, an imaging apparatus is provided. The 



[0006] In a further aspect, a computer readable medium is provided 
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determine a attenuation difference X(k) in accordance with A(k) = r(k) - , 
calculate an amount of attenuation to add r|(k) in accordance with 

rj(k) = y Rf (k) - x t (k)pi (k) - Rf (k) arcsinj^^ j + 1 Rf (k) - x r (k)p r (k) - Rf arcsin 

where pi(k), and p r (k), are the magnitude of a left and a right projection boundary 
samples averaged over multiple views, respectively, and xi(k), x r (k), Ri(k), and R r (k) 
are a location and radius of a left and right fitted cylinders, respectfully, compare r|(k) 

Tj( k} Ll 

to X(k) by calculating a ratio s{k) = — — where jx w is an attenuation coefficient of 

A(k) 

water, compare e(k) to a threshold q, and use at least one of r|(k) and X(k) to correct 
an image when e(k) is not greater than q. The program is further configured to 
instruct the computer, when e(k) is greater than q, to calculate a r| n (k) based on data 
regarding a k t -n view and a k2+n view, wherein n is an integer, and correct an image 
using the Ti n (k). 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] Figure 1 is a pictorial view of a CT imaging system 

embodiment. 

[0008] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

[0009] Figure 3 illustrates a projection truncation occurring only in a 
portion of the projection angles. 

[0010] Figure 4 illustrates a reconstructed image in which no 
truncation correction was performed. 

[0011] Figure 5 illustrates an image with truncation correction as 
described in attorney docket number 129993. 
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[0012] Figure 6 is a graph of a total amount of measured attenuation 
of each projection plotted as a function of projection view. 

[0013] Figure 7 graphically describes a fitting process. 

[0014] Figure 8 illustrates a plurality of images including an image 
where only a single view was used to estimate the magnitude and the slope of the left 
truncated projections, pi, S], and a single view was used to estimate the magnitude and 
the slope of the right truncated projections p r , and s r . 

[0015] Figure 9 graphically illustrates how missing projection 
samples are calculated based on the estimated cylinders. 

[0016] Figure 10 shows the same phantom scan shown in Figures 4 
and 5 corrected with the herein described methods and apparatus. 

DETAILED DESCRIPTION OF THE INVENTION 

[0017] There is herein provided truncation compensation methods 
and apparatus for Extended Field-of-View in Computed Tomography Systems. The 
apparatus and methods are illustrated with reference to the figures wherein similar 
numbers indicate the same elements in all figures. Such figures are intended to be 
illustrative rather than limiting and are included herewith to facilitate explanation of 
an exemplary embodiment of the apparatus and methods of the invention. 

[0018] In some known CT imaging system configurations, a radiation 
source projects a fan- shaped beam which is collimated to lie within an X-Y plane of a 
Cartesian coordinate system and generally referred to as an "imaging plane". The 
radiation beam passes through an object being imaged, such as a patient. The beam, 
after being attenuated by the object, impinges upon an array of radiation detectors. 
The intensity of the attenuated radiation beam received at the detector array is 
dependent upon the attenuation of a radiation beam by the object. Each detector 
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element of the array produces a separate electrical signal that is a measurement of the 
beam attenuation at the detector location. The attenuation measurements from all the 

detectors are acquired separately to produce a transmission profile. 

/ 

[0019] In third generation CT systems, the radiation source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged such that an angle at which the radiation beam intersects the object 
constantly changes. A group of radiation attenuation measurements, i.e., projection 
data, from the detector array at one gantry angle is referred to as a "view". A "scan" of 
the object includes a set of views made at different gantry angles, or view angles, 
during one revolution of the radiation source and detector. 

[0020] In an axial scan, the projection data is processed to 
reconstruct an image that corresponds to a two dimensional slice taken through the 
object. One method for reconstructing an image from a set of projection data is 
referred to in the art as the filtered back projection technique. This process converts 
the attenuation measurements from a scan into integers called "CT numbers" or 
"Hounsfield units", which are used to control the brightness of a corresponding pixel 
on a display device. 

[0021] To reduce the total scan time, a "helical" scan may be 
performed. To perform a "helical" scan, the patient is moved while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from a 
fan beam helical scan. The helix mapped out by the fan beam yields projection data 
from which images in each prescribed slice may be reconstructed. 

[0022] As used herein, an element or step recited in the singular and 
proceeded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
to "one embodiment" of the present invention are not intended to be interpreted as 
excluding the existence of additional embodiments that also incorporate the recited 
features. 
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[0023] Also as used herein, the phrase "reconstructing an image" is 
not intended to exclude embodiments of the present invention in which data 
representing an image is generated but a viewable image is not. Therefore, as used 
herein the term, "image," broadly refers to both viewable images and data 
representing a viewable image. However, many embodiments generate (or are 
configured to generate) at least one viewable image. 

[0024] Figure 1 is a pictorial view of a CT imaging system 10. 
Figure 2 is a block schematic diagram of system 10 illustrated in Figure 1. In the 
exemplary embodiment, a computed tomography (CT) imaging system 10, is shown 
as including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has a radiation source 14 that projects a cone beam 16 of X-rays toward a 
detector array 18 on the opposite side of gantry 12. 

[0025] Detector array 18 is formed by a plurality of detector rows 
(not shown) including a plurality of detector elements 20 which together sense the 
projected X-ray beams that pass through an object, such as a medical patient 22. Each 
detector element 20 produces an electrical signal that represents the intensity of an 
impinging radiation beam and hence the attenuation of the beam as it passes through 
object or patient 22. An imaging system 10 having a multislice detector 18 is capable 
of providing a plurality of images representative of a volume of object 22. Each 
image of the plurality of images corresponds to a separate "slice" of the volume. The 
"thickness" or aperture of the slice is dependent upon the thickness of the detector 
rows. 

[0026] During a scan to acquire radiation projection data, gantry 12 
and the components mounted thereon rotate about a center of rotation 24. Figure 2 
shows only a single row of detector elements 20 (i.e., a detector row). However, 
multislice detector array 18 includes a plurality of parallel detector rows of detector 
elements 20 such that projection data corresponding to a plurality of quasi-parallel or 
parallel slices can be acquired simultaneously during a scan. 
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[0027] Rotation of gantry 12 and the operation of radiation source 14 
are governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes a radiation controller 28 that provides power and timing signals to radiation 
source 14 and a gantry motor controller 30 that controls the rotational speed and 
position of gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 
samples analog data from detector elements 20 and converts the data to digital signals 
for subsequent processing. An image reconstructor 34 receives sampled and digitized 
radiation data from DAS 32 and performs high-speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

[0028] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, 
radiation controller 28, and gantry motor controller 30. In addition, computer 36 
operates a table motor controller 44 which controls a motorized table 46 to position 
patient 22 in gantry 12. Particularly, table 46 moves portions of patient 22 through 
gantry opening 48. 

[0029] In one embodiment, computer 36 includes a device 50, for 
example, a floppy disk drive or CD-ROM drive, for reading instructions and/or data 
from a computer-readable medium 52, such as a floppy disk or CD-ROM. In another 
embodiment, computer 36 executes instructions stored in firmware (not shown). 
Generally, a processor in at least one of DAS 32, reconstructor 34, and computer 36 
shown in Figure 2 is programmed to execute the processes described below. Of 
course, the method is not limited to practice in CT system 10 and can be utilized in 
connection with many other types and variations of imaging systems. In one 
embodiment, Computer 36 is programmed to perform functions described herein, 
accordingly, as used herein, the term computer is not limited to just those integrated 
circuits referred to in the art as computers, but broadly refers to computers, processors, 
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microcontrollers, microcomputers, programmable logic controllers, application 
specific integrated circuits, and other programmable circuits. Although the herein 
described methods are described in a medical setting, it is contemplated that the 
benefits of the invention accrue to non-medical imaging systems such as those 
systems typically employed in an industrial setting or a transportation setting, such as, 
for example, but not limited to, a baggage scanning CT system for an airport or other 
transportation center. 

[0030] Some known methods have been published that address the 
artifact reduction but not the imaging of a portion of patient 22 that is outside a field 
of view (FOV). However, it is desirable to image the portion of the patient that 
extends beyond the FOV. This is useful in many fields including Oncology, Spin 
Angiography, Fused imaging systems, and in Economy CT Scanners. The current 
hardware of at least one known multi-slice CT scanner limits the reconstruction field 
of view (FOV) to about 50 centimeters (cm.). Although this is sufficient for most 
clinical applications, it is desirable to expand the FOV to image objects outside this 
FOV. This may have particular advantages for applications such as oncology or 
CT/PET. For oncology apphcations, a larger FOV is desired. This is mainly due to 
the fact that for radiation treatment planning, the limbs of the patient are often 
positioned outside the scan FOV for better tumor positioning. The known CT 
reconstruction algorithms ignore the truncated projections and produce images with 
severe artifacts. These artifacts may affect an accurate estimation of the attenuation 
path for treatment planning. Herein described is an algorithmic approach to increase 
the reconstruction FOV beyond the FOV limited by the detector hardware. This 
correction algorithm can be applied to various reconstruction algorithms including, 
but not limited to, Full scans, Halfscan/segment, helical scan, and cardiac sector based 
algorithms. Additionally, system 10 is configured to employ the herein described 
algorithms. 

[0031] Figure 3 illustrates a projection truncation occurring only in a 
portion of the projection angles. Under some scanning conditions, portions of patient 
22 extend beyond the region measured by detector 18, which may lead to image 
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artifacts and an incomplete representation of the imaged object. The X-ray tube and 
detector 18 are rigidly held together on a frame which rotates around patient 22. 
During the course of rotation, measurements are continuously made within a "fully 
sampled field of view" 60. The attenuation of x-rays traversing any areas of object 22 
positioned outside fully sampled field of view 60 is measured at a limited range of 
rotation angles, and this region is referred to as the "partially sampled field of view" 
region. In other words, portions that are within fully sampled field of view 60 are 
positioned within fan 16 such that measurements are obtainable at all gantry angles, 
and the data collected is defined as fully sampled field of view data. Some portions, 
however, are within fan 16 at some angles but are outside fan 16 at other angles, and 
the data collected regarding these portions is defined as partially sampled field of view 
data. In this example, the projection taken at a 3 o'clock position is free of truncation 
and the projection taken at a 12 o'clock is severely truncated. Therefore, one can rely 
on the un-truncated projections (i.e., positions around 3 o'clock in Figure 3, for 
example) to estimate the amount of truncation for the truncated views (e.g., positions 
around 12 o'clock for the example in Figure 3). One early step in the correction 
process is to perform software fan-beam to parallel-beam rebinning on the pre- 
processed projections. In one embodiment, that early step is the first step. This 
process is well known in the art and does not require special data collection. Once the 
rebinning is completed, the projections are integrated over all detector channels to 
obtain a total attenuation curve as described in more detail below. 

[0032] Figure 4 illustrates a reconstructed image 80 in which no 
truncation correction was performed. As is easily observed, image 80 contains 
truncation artifacts 82. 

[0033] Attorney docket number 129993 describes an algorithm that 
enables truncation compensation for a majority of situations. Under certain 
conditions, however, the compensation algorithm described in attorney docket number 
129993 does not produce optimal results. Figure 5 illustrates an image 84 with 
truncation correction as described in attorney docket number 129993. Although 
image 84 illustrates much improvement over image 80, image 84 has distortion in a 
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corrected area 86 (restored arm). In this example, a chest phantom was positioned so 
that an outermost region is 29cm from iso-center (at the edge of a 58cm FOV). 
Although the previous correction removes the artifacts inside the 50cm FOV, the 
shape and intensity of the restored arm 86 is somewhat distorted. 

[0034] It is believed that the distortion is caused mainly by an 
inaccuracy in the estimation of the missing attenuation. For illustration, a total 
amount of measured attenuation 90 of each projection is plotted as a function of 
projection view in Figure 6. For clarity, the graph is plotted on an expanded scale. A 
plurality of dips 92 in measured attenuation 90 correspond to truncated views. It is 
clear from the graph that the summation of the measured projection measurement for 
each projection is not a constant over projection angles for a plurality of non-truncated 
projections (areas of measured attenuation 90 which are approximately flat). This is 
due to a non-ideal calibration and the fan-parallel rebinning process. Although the 
variation is small (less than 2%), it could represent a significant portion of the total 
truncated projection samples. In other words, as compared to the estimated missing 
projection, the variation percentage is much larger. Therefore, using additional 
information for the truncation estimation than used in the compensation described in 
attorney docket number 129993 will result in less or no distortion. 

[0035] During a plurality of phantom experiments, it was found that 
the magnitude and slope at the projection boundary alone are sometimes insufficient 
for the estimation of the size and location of a cylindrical object to be fitted in to the 
projection because sometimes the slope of the boundary sample(s) is biased by an 
internal structure of the scanned object. For example, if an air pocket happens to be 
located tangential to a truncated projection boundary sample, significant discontinuity 
(in projection angle) in the extrapolated projection samples can result. This could lead 
to shading and streaking artifacts. 

[0036] Another complication occurs when the projection samples are 
noisy (due to limited x-ray photon statistics). Under this condition, both the 
magnitude and slope of the projection boundary may become unreliable. High 
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frequency variations (from view to view) may result in severe streaks and highly 
distorted object shapes. Clinically, this can occur when a large patient is scanned or 
when a low-flux acquisition protocol is used. The herein described methods and 
apparatus at least partially overcome these above described difficulties. 

[0037] More specifically, projections are declared either truncated or 
non-truncated (i.e., un-truncated). To declare (i.e., determine) a projection as 
truncated, a condition is imposed that either side of the projection boundary samples 
(averaged over a few channels) has to be less than a pre-determined threshold in 
addition to the previously imposed conditions identified in docket number 129993. 
The average of plurality of samples should be an average of plurality of edge samples. 
That is, if the projection contains N samples indexed from 0 to N-l, the average of the 
left m projection samples from 0 to m-1, and the average of the right m samples from 
N-m to N-l are used. In addition to the above test, the test in attorney docket 129993 
is still used. That is, the sum of all N samples is calculated for all views. The 
maximum value over all projection views is calculated, and when the sum of a 
particular view is less than a pre-determined percentage of the maximum value, this 
view is declared truncated. The threshold value is non-zero to account for the non- 
ideal calibrations. Mathematically, a projection is declared truncated if the following 
condition exists 



m m m 



Pi =— yp(i,*)>/,w p r = — S^p(N-i,k)>t 



m m 
i=l i=l 

[0038] where N is the number of detector channels, k is the 
projection view index, m is the number of samples, and t is the value of the threshold. 
In an exemplary embodiment, t equal to about 0.42 has empirically shown to provide a 
good performance. In other embodiments, t is between about 0.25 and 0.58, t is 
between about 0.33 and 0.5, and/or t is between about 0.375 and 0.46. 

[0039] Using this above described determination method, the non- 
truncated views are identified as shown with double arrows in Figure 6. To estimate 
the amount of missing projections, the total attenuations of neighboring un-truncated 
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projections are used for the estimation. In specific, the total attenuation of the 
truncated view is obtained by linearly interpolating the two un-truncated projections 
on both sides of the truncation view region. For the example shown in Figure 6, the 
first set of truncated views is located between projection view number 290 (106° in 
terms of projection angle) and 445 (163° in terms of projection angle). The total 
attenuation obtained for view number 289 and 446 is used to estimate the total 
attenuation of the truncated views. The estimation is performed by linearly 
interpolating the two un-truncated samples, as shown by dotted lines in Figure 6. 
Mathematically, the estimated total attenuation, x(k), is then 

r<*) = j^«*l) + -p^£(*2) ( 2 ) 

[0040] where ki and k2 are the view locations of the un-truncated 
views adjacent to the truncation region, ^(k) is the total attenuation obtained by 
summing up the projection measurement for view k. Mathematically, it is described 
by the following equation 

N 

m = ^P(Uk) (3) 

[0041] To improve the robustness of the algorithm, one can use the 
average of several projection views for the un-truncated projections. For example, 
one could use the average total attenuation of views 285 to 289 instead of the view 
289 alone, and ki then refers to the average values for views 285 to 289. In an 
exemplary embodiment, the average is taken over five views, but in other 
embodiments, the average is taken over a number of views other than five (e.g., 3-7 
views, and 2-8 views). Additionally, k 2 can refer to a single view or an average of 
several view independent of ki . For example, and referring the to example shown in 
Figure 6, let ki refer to the average values of views 285 to 289, then, in one 
embodiment, k2 refers to the average values of views 446-450, and, in another 
embodiment, k2 refers to 446. Additionally, ki is independent of k 2 , and ki can refer 
to a single view or multiple views regardless of k 2 . 
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[0042] An additional amount of projection needed to make up for the 
truncation, X(k), is then the difference between estimated total attenuation and the 
measured total attenuation 

= *(*)-£(*) (4) 
[0043] Next, a size and a location of a water cylindrical object that 
can be fitted at the boundaries of each truncated projection are estimated to come up 
with an estimated projection. The location and the radius of the fitted cylinder is 
calculated based on the magnitude and the slope of the truncated projection boundary 
samples averaged over m projection views 

. and R m . Jfl!ffi + 2 W (5) 

^ (t) = Sr(k)Pr(k) ^ ^ ^ {k) = + ^ 2 (Jfc) (6) 

[0044] where pi, s\ 9 p r , and s r are the magnitude and slopes of the left 
and right projection boundary samples averaged over multiple views, respectively, and 
xi, x r , Rj, and R r are the location and radius of the left and right fitted cylinders, 
respectfully, and ja w is the attenuation coefficient of water. The fitting process is 
described graphically in Figure 7. 

[0045] Previously, it was indicated that under x-ray photon starvation 
conditions, the parameter estimation process might become unreliable because the 
estimated parameters may fluctuate significantly from view to view. Figure 8 
illustrates a plurality of images 100 including an image 102 where only a single view 
was used to estimate the magnitude and the slope of the left truncated projections, pi, 
Si, and a single view was used to estimate the magnitude and the slope of the right 
truncated projections p r , and Sr. Image 102 clearly illustrates streaking artifacts caused 
by parameter fluctuation. To overcome this possible shortcoming, the use of 
projection samples averaged over several views is employed, in one embodiment, to 
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estimate the magnitude and the slope of the truncated projections, pi, si, p r > and s r . In 
an exemplary embodiment, projection samples of seven views are used for the slope 
and magnitude estimation. Plurality of images 100 also includes an averaged image 
104 is which seven views were used for the slope and magnitude estimation. In 
another embodiment, between 5 and 9 views inclusive are used for the slope and 
magnitude estimation. The averaged values are then used for the radius and location 
determination of the fitted cylinders. Based on the fitted cylinders, the total amount of 
attenuation of the added cylinders for projection k, r|(k), (the portion of the cylinder 
that extends beyond the original projection) is calculated based on the following 
formula 



estimated amount of additional attenuation, X(k), described in Eq. (4). Ideally, these 
two quantities should be equal. In reality, however, since the missing portion of the 
object is rarely a perfect cylinder, there will always be some differences between the 
two quantities. Therefore, the ratio of the two quantities is calculated, e(k) 



potentially indicates a problem with the fitted cylindrical parameters. This could be 
caused by an irregular shape of the scanned object as described previously. Therefore, 
when s(k)>q, one may not want to use the estimated parameters for projection 
compensation, q is a pre-determined threshold and is set to 2 in an exemplary 
embodiment. In another exemplary embodiment, q is between 1.5 and 2.5. In another 
embodiment, q is between 1.75 and 2.25. In another exemplary embodiment, q is 
between 1.9 and 2.1. 





[0047] When the quantity, s(k), is significantly different than unity, it 



-14- 



Docket No. 134692 



Express Mail No. EV 339984887 US 



[0048] For the projections in which e(k)>q, one may want to use 
interpolated values of the cylindrical parameters of two projections with a good 
parameter fit (i.e., e(k)<q). That is, if e(k) values for views ki to k2 exceed the 
threshold, one may want to use the fitted cylinder parameters of view ki-n and k2+n to 
estimate the cylindrical parameters for these views, where n is a parameter, and n=l 
indicates the first "good fit" views. In an exemplary embodiment, n=5 was selected 
and linear interpolation has been shown to be sufficient. That is, xj(k), x r (k), Ri(k), 
and R r (k) are linearly interpolated based on xj(kl-n), xj(k2+n), x r (kl-n), x r (k2+n), 
Ri(kl-n), Rj(k2+n), R r (kl-n) and R r (k2+n), and r| n (k) is calculated according to Eq. 7 
using the interpolated values. The ratio, s(k), is also linearly interpolated from the 
same set of views. Of course, other interpolation method, such as higher order 
interpolations, can be used. In addition, additional low-pass filtering along the view, 
k, can be applied to the estimated parameters to further improve the continuity of the 
parameters. In other embodiments, n is between 2 and 8, and/or n is between 3 and 7. 

[0049] Finally, the missing projection samples are calculated based 
on the estimated cylinders. Figure 9 graphically illustrates this process. Since the 
total attenuation of these cylinders is unlikely to match the total missing attenuation 
estimated using the consistency condition (Eq. 4), additional adjustment can be used. 
For example, the x-axis (detector channel index i) is scaled by the ratio so that the 
initial estimated projection (shown by a dotted line in Figure 9) is further expanded 
(shown by a solid thick line in Figure 9) or compressed in x. This can be 
accomplished by changing the indexing distance of detector channel i or by 
interpolating the initially calculated projections. 

[0050] With the herein described methods and apparatus, the 
performance and robustness of the methods and apparatus described in attorney 
docket number 129993 are significantly improved. Figure 10 shows the same 
phantom scan shown in Figures 4 and 5 corrected with the herein described methods 
and apparatus. Improvement in the accuracy of the shape and density of an object 106 
outside the scan FOV is obvious, which is one technical effect or contribution for the 
disclosed invention. Of course, other technical effects exist also. 
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[0051] Although the above described system and methods uses only 
the conservation of total attenuation, the magnitude, and the slope of the boundary 
samples to estimate the missing projection distribution, additional information can 
also be used for the estimation. For example, one could use the Helgason-Ludwig 
Condition (HL condition) for tomography to further refine the above techniques. 
Additionally, different thresholds can be placed to ensure that the algorithm functions 
properly under erroneous measurement conditions. For example, one could set the 
upper and lower limits on the stretch ratio described in Figure 9 to prevent the 
condition of increased error due to unreliable measurement. In addition, the slope 
calculation of Si and Sr can be set such that it falls within a reasonable range. If the 
characteristic of the material of the scanned object is known to be significantly 
different from water, one can also use the attenuation coefficients of the known 
material (instead of water) to perform size and location calculations shown in Eqs. (3) 
and (4). 

[0052] Because the interpolated data may not have the same image 
quality as data within the fully sampled FOV it may be useful to tag the image where 
the FOV becomes extrapolated. In one embodiment, a delineation is provided in a 
reconstructed image between areas representative of the fully sampled field of view 
data and the partially sampled field of view data. This could also be done with a color 
code or shift in the CT Number. Because the tag may impact the ability to view the 
image data an easy way is provided to turn on and off the tagging. And therefore, a 
user of system 10 is allowed to turn on or off the tagging. 

[0053] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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